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Abstract 



We study the level-statistics of a disordered system undergoing the Ander- 
son type metal- insulator transition. The disordered Hamiltonian is a sparse 
random matrix in the site representation and the statistics is obtained by 
taking an ensemble of such matrices. It is shown that the transition of lev- 
els due to change of various parameters e.g. disorder, system size, hopping 
rate can be mapped to the perturbation driven evolution of the eigenvalues 
of an ensemble subjected to Wigner-Dyson type perturbation with an initial 
state given by a Poisson ensemble; the mapping is then used to obtain desired 
level-correlations . 

PACS numbers: 05.45+b, 03.65 sq, 05.40+j 
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. The degree of disorder and dimensionality of a disordered system plays a significant 
role in the metal- insulator (MI) transition [|TJ. The level statistics during the transition is 
governed by a parameter which is a function of both, disorder as well as dimensionality [0 . 
The purpose of this paper is to identify the parameter and study the effect of its variation 
on the spectral fluctuations. 

The statistical behaviour of the levels of a disordered Hamiltonian depends on their 
number r within the energy ranges of interest and the dimensionless conductance g. The 
parameter g is the ratio of two energy scales, namely, Thouless energy E c and mean level 
spacing A. The energy scale E c = hD/L 2 , with D as diffusion coefficient and L system size, 
corresponds to ergodic time-scales (ergodic classical limit) above which the wavefunction 
has fully diffused in the entire volume of the system and therefore does not feel the space 
dimensionality d. The eigenfunctions associated with levels within energy ranges smaller 
than E c are extended structureless objects leaving no signature on the level correlations. 
The number r < g, therefore, results in a universality of level correlations which can be well- 
modeled by the invariant ensembles of random matrix theory (RMT) ||. For r > g, the time 
scales associated with the dynamics are such that the wavefunction is not yet fully diffused 
in the whole volume. The eigenfunctions in the corresponding energy-range are localized, 
either partially or fully depending on the dimensionality, and their spatial correlations can 
significantly affect the level correlations. 

The metallic phase of a disordered system corresponds to g — > oo in the thermodynamic 
limit L — > oo. The states, therefore, in any large but finite interval are always extended 
with a Wigner-Dyson (WD) level-statistics, characterstic of strong level- repulsion ||. In the 
insulator phase g — > in the thermodynamic limit, implying no overlap between eigenfunc- 
tions. The corresponding level statistics can be modeled by the statistics of independent 
random numbers, that is, Poisson statistics. A change in the degree of disorder results in the 
variation of g (for fixed r) which can induce a transformation of extended states into local- 
ized states thereby causing a transition of spectral fluctuations. However the MI transition 
is very sensitive to the dimensions of the system. For d — 1,2 and, in thermodynamic limit, 
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the system always remains an insulator even for very weak disorder. For d > 2, the increase 
of disorder results in decreasing g, leading to an extended to localized state transition at a 
critical value of disorder. In the region near the critical point of the transition, known as 
Anderson transition, the eigenstates are neither localized nor extended but form self-similar 
multifractal structures with strong fluctuations at all length scales ||. This results in a 
new type of universal level-statistics, known as critical spectral statistics (CSS), which is 
neither Poissonian nor Wigner-Dyson but is a hybrid of both For finite system size 

also, a crossover from Poisson to WD statistics can be seen among levels as a function of the 
parameter (/L with ( as the correlation length. However the detailed analytical information 
about various level correlations in the critical regime for inifinte systems as well as for the 
intermediate stages, during transition in finite systems, is still missing. Our study attempts 
to bridge this gap. 

The level statistics for various transition stages is significantly affected by the nature of 
the localization of the eigenfunctions and can suitably be modeled by an ensemble containing 
eigenvalue-eigenfunction correlations (therefore basis-dependent). Recently a new technique 
has been developed to deal with such cases where the eigenvalue distribution for an ensemble, 
non-invariant under a change of basis, is mapped to the one for an ensemble invariant under 
change of basis and subjected to a random perturbation ||. The mapping is achieved by 
identifying a basis-dependent parameter in the former to the perturbation (e.g symmetry- 
breaking) parameter in the latter. The distribution for the invariant case can subsequently be 
mapped to the non- stationary states of the particles governed by Calogero-Sutherland-Moser 
(CSM) Hamiltonian ||. The known particle correlations can therefore be used to determine 
level- correlations of the non-invariant ensembles. The present study uses the technique to 
analyze the spectral properties of the levels undergoing MI transition by taking an ensemble 
of Anderson Hamiltonains in the site representation. As expected, the basis dependent 
parameter, governing the transition in this case, is a function of disorder, dimensionality 
and hopping rate and is related to g. 

The Anderson model for a disordered system is described by a d-dimensional disordered 
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lattice, of size L, with a Hamiltonian H = J2n e nat a n + J2n^mt( a i a m + dnO-m) m tight- 
binding approximation [I]. Here a+ and a n are the creation and annihilation operators of 
an electron at a site n in a lattice and m refers to the nearest neighbors of the site n. The 
site energies e„, measured in units of the overlap integral between adjacent sites correspond 
to the random potential. The hopping is generally assumed to connect only the z nearest- 
neighbors with amplitude t so that the electron kinetic energy spread or bandwidth is zt. In 
the configration space representation the Hamiltonian H turns out to be a sparse matrix of 
size N = L d with diagonal matrix elements as the site-energies e«. The level-statistics can 
therefore be studied by analyzing the properties of an ensemble of (i) sparse real symmetric 
matrics in presence of a time-reversal symmetry and (ii) sparse complex Hermitian matrice 
in absence of a time-reversal symmetry. 

We first consider the case of MI transition brought about by decreasing the diagonal 
disorder. In this case, site-energies e« are taken to be independent random variables with 
probability-density p(ej). In the original Anderson model |I] p(e) was taken to be a constant 
W~ l between —W/2 to W/2. Various physical arguments and approximations used in this 
case led to conclusion that all the states are localized for W > iKtln(^) with K as a 
function of z and d. However, as well-known now, MI transition does not depend on the 
nature of p(e) and latter can also be chosen as Gaussian; the type of p(e) affects only 
the critical point of the transition [0. The ensemble measure p(H), for any intermediate 
state of MI transition brought about by diagonal disorder, can therefore be described by 
p(H, y, b) oc exp(- Y? s=1 Efc at kk . s (H kk . s - b kk . s f) Y[ k <i-s=i,2^{Hu;s ~ hi;s) with subscript "s" 
referring to real and imaginary parts and the Kronecker delta function implying the non- 
randomness of off-diagonal matrix element (H k i = b k i-i + ib k i-2)- However p(H) can be 
expressed in a more general form: 
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p(H, y, b) = Cexp(- a ki;s(H k i ;s - b H . s ) 2 ) (1) 

s=l k<l 

with C = T\ k <i Ils=i y^§^ as the normalization constant, y as the set of the coefficients 
Uki,s = atki;s9ki = iJirl > and b as tne set of all b H . s with g kl = 1 + 5 ki . As obvious, in limit 
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otki-,1, &ki;2 — * oo, eq.(l) corresponds to the non-random nature of the off-diagonal part of 
H. However note, for later reference, that a same mean value for all the matrix elements, 
that is bki- s = e for all k, I (irrespective of a.ki values) implies just a shift of the origin in the 
matrix space and does not affect the eigenvalue distribution. 

The MI transition is brought about by a competitive variation of the disorder and hopping 
rate which can be mimicked by a change of the distribution parameters of p(H). The changed 
eigenvalue distribution as a result can be obtained by integrating p(H) over associated 
eigenvector space. Let P(p, y, b) be the probability of finding eigenvalues Aj of H between 
Pi and pi + dpi at a given y and b, it can be expressed as P(p,y,b) = JY[iLi$(Pi — 
\i)p(H, y, b)dH. As discussed in ref. [[|, a particular combination of the evolution of P with 
respect to various distribution parameters, namely, DP = 2 J2 S J2k<i Vki;s{l — Vm-,.^ '''' 



dVk 



lUs Hk<i bki sj^ — , with 7 as an arbitrary parameter, results in the diffusion of eigenvalues 
with a finite drift due to their mutual repulsion. 
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with (3 depending on the symmetry class of the ensemble ((3 = 1 for real symmetric matrices 
and (3 = 2 for complex Hermitian ones). The multi-parametric evolution of P can further 
be reduced in terms of the evolution with respect to a single parameter Y, that is, DP = |^ 



where Y is a function of various y^i-s and bki- s 
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and can therefore be referred as the 'complexity parameter'. Here C is an arbitrary 
constant. The eq.(2) can now be rewritten as follows, 
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For yki-s — > 7 and bki- s — * e (for almost all k, I), DP — > and the evolution reaches the 
steady state given by a WD distribution P(p) = ]lj<j I A** — Pj\^ e ~^ (Note, in this 

limit, p(H) oc e -7TrH 2 -eTr// ^jp^ indeed has a WD distribution and therefore agrees with 



the solution of eq.(4) in limit DP — > 0). The eq.(4) is similar to the equation governing a 
perturbation driven transition of the eigenvalues, of a Hamiltonian H = H + tV, between 
various universality classes of WD ensembles. Here Hq is an initial ensemble with pertur- 
bation V causing a transition to another (or same) universality class. The transition to 
equilibrium, with perturbation strength r as the evolution parameter, is rapid, discontinu- 
ous for infinite dimensions of matrices (analogous to Anderson transition). But for small-r 
and large N, a smooth transition can be seen in terms of a rescaled parameter A which 
measures locally the mean-square perturbation matrix element in units of the average level 
spacing [0. The intermediate states of this transition have been well-studied in past and 
many results for their spectral fluctuations are already known By replacing the per- 

turbation parameter by complexity parameter, these results can directly be used for the 
corresponding measures for the MI transition. Note this is eqivalent to say that if the initial 
Anderson Hamiltonian is described by Hq, the changing complexity (the combined effect 
of changing disorder and hopping rate) can be treated as a perturbation V of strength Y, 
taken from WD ensemble, with Hamiltonian H of the system represented by H = H + YV. 

The n th order level density correlation R n is defined as R n =< v([/,i, F)..z/(/z n , Y) > with 
z/(/i, Y) = A -1 J2i — Hi) as the density of eigenvalues and < .. > implying the ensemble 
average. The R n can also be expressed in terms of P: R n = n^z^i I P(fa ^)d/i n +i..d/xjv- A 
solution of eq.(4), averaged over initial conditions, can therefore be used to obtain various 
correlations for a given set of distribution parameters. In practice, however, the involved 
integrals could only be solved for the cases when H is complex Hermitian in nature a 
direct integration of F-P equation (4) leads to the BBGKY hierarchic relations among the 
unfolded correlators Rnin, .., r n ; A) = LimN -> oo; r^.yY'r^Iy) with r = f R i(^ Y ) dY 
and A = (y — loh/A 2 (A as the mean level spacing = R^ 1 = (N < v >)~ x ) g,|], 

dA y dr) drj \rj - r k J y drj J-oo r s - r k 

Here note a perturbation of strength A mixes levels in a spectrum span of v^A only, due 
to level repulsion 0. This results in a sensitivity of the correlations within a given energy 
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range r to the ratio r/\fX. In the case of MI transition, therefore, A can be identified with 
9 2 - 

For n = 2 and small values of r, the integral term in eq.(5) makes a negligible contribution 
thus leading to following approximated closed form equation for R 2 

^2 = 2 ^_ 2 ^ ( 6) 
<9A dr 2 dr r 

Similarly for large-r behaviour, R 2 can be obtained from the following relation 

/ ,x ^ / \ n „ k f°° , i?2(^ - s; 0) - i? 2 (r - s; oo) 
R^K) = R 2 (r^ )+ 20Kj_J s >- (7) 

The hierarchic equation can then be used to obtain an approximate form of the higher order 
correlations. For example, the approximate information about R3 can be extracted by a 
substitution of large and small r behaviour of R 2 in eq.(5) with n — 2. 

For studies of the level-correlations during MI transition, the system can initially be 
considered in an insulator regime where all the eigenvectors become localized on individual 
sites of the lattice (strong disorder limit). This results in a diagonal form of the matrix 
H with the eigenvalues independent from each other. The insulator limit can therefore 
be modeled by ensemble (1) with a^i — > 00 for k 7^ /, = aj (for all ^-values) and 
Hi ~^ e (f° r ah k, I), giving, Y = ^-jyln ^""^ — ^lne + C. (Here the origin of the energy-axis 
is chosen at e so as to avoid the undefined nature of ln6 w near b ki — > 0). The decrease 
of the diagonal disorder, that is, an increase of a k k from a to some finite values (while 
aw, k 7^ I, remains infinite throughout the transition) will ultimately lead to metal regime 
with fully delocalized wavef unctions. A change in the hopping rate of the connected sites 
from e to e + t (pki — > e + t for only those k, I values which are connected, bki = e for 
others) will also affect the transition. As mentioned above, the eigenvalue distribution of 
H in the metal regime can be well-modeled by the WD ensemble; let it be described by all 
ttfcz — ► «m(> cco). Thus for the study of transition in this case we can choose 7 = 2a m- For 
any intermediate state of the transition with hopping rate t and disorder coefficient a, the 
complexity parameter Y = ^2 [f m p|z^ - f + e ) - 2N2 ~ K \ne] + C. Here K = kN is 



the total number of the sites connected by hopping and depends on the dimensionality d of 
the system (as iV = L d ). The transition parameter can now be given as follows, with the 
mean level spacing A = (< v > N)' 1 : A = 2a M (Y - Y )/A 2 = a M < v > 2 Na(a,t)/2 with 

a(a,t) = \ n a \ a °~~ aM \ !^i n t±i _ As obvious from the above, the transition is governed 

by relative values of the disorder and the hopping. Here A — > leads to a fully localized 
regime which, corresponds to the condition a(a,t) < ■h. Similarly the condition a(a,t) > 
corresponds to A — > oo for large A" and therefore extended states. For a(a,t) = (/i as 
a ^-independent function), A = % < v > 2 = A* is independent of A" which therefore 
remains same for infinite systems (< v > being independent of N). This A- value corresponds 
to the critical point of the transition in infinite systems, the statistics being Poisson and WD 
below and above A*, respectively. The condition for the critical region or mobility edge can 
therefore be given as hi~ + "~ A " Q ~ ^j^ 11 (^ + 1) (taking e = 1 without any loss of generality) . 

As ][a ~ a ° ]{ << 1 even for large a- values, the condition is always satisfied if — > 0. This 

explains the localization of all the states in infinitely long wires (or strictly 1-d systems where 
z is very small) even for very weak disorder. With increasing dimensionality d, connectivity 
K of the lattice and thereby the possibility of A >> and the delocalized states increases. 
As reflected by its form, A also plays the role of the finite size scaling parameter. Note that 
the A-dependence of various correlators indicates an interaction of levels different from that 
of Poisson or WD if A remains finite in limit A^ —>■ oo. This implies the presence of a new 
type of universal statistics near the critical point A* of the transition. For finite systems, 
however, there is a continuous family of intermediate statistics between Poisson and WD, 
described by A and the symmetry of the Hamiltonian. 

The fluctuation measures for various stages of the MI transition can be obtained from 
those of the perturbation driven Poisson — > WD ensemble transition, by replacing perturba- 
tion parameter by Y. The two level density correlator i?2( r ; A) for the Anderson transition 
in finite systems, in presence of a magnetic field, can therefore be given by the R2 for Poisson 
— > GUE transition driven by the perturbation P,pH|]: 



S 



4 roo rl 

R 2 (r; A) — R 2 (r; oo) = — dx dz cos(27irx) exp — 87r 2 Aa;(l + x + 2z\fx) 

n Jo J-l L 



z 2 ){l + 2z^x~) 



1 + X + 2Z\fx 

where R 2 (r, oo) = 1 — s ™ 2 ^ r ' ) and -R2(f, 0) = 1 corresponding to metal and insulator regime 
respectively. A substitution of A = A* in eq.(8) will thus give the two level correlation for 
critical regime of Anderson transition. 

An important characterstic of CSS is the level compressibility \ which is basically a 
measure of the level repulsion; x = in the metallic phase and x — 1 i n the insulator phase 
and takes an intermediate value between and 1 at the hybrid phase near the critical point. 
The eq.(4.15) of can be used to obtain x(A) « 1 - J^° 00 Y 2 (r; A)dr with Y 2 = 1 - R 2l 



X(A) = l + J oo du^[I 1 (z)-^8z/7iI 2 (z)}exp[-(A/2)u 2 -n^A/2u] with z = V2 * ttA 2 m 3 and 
/„ as the n th Bessel function. (Note eq.(8) can not directly be used to calculate x due to 
technical difficulties). It can be checked that x — > 1 for A — ► (insulator limit) and x ~ > 
for A — > oo (metal limit). As obvious from the smooth nature of the functional form of x-, 
any non-zero finite A-value will correspond to a value of % intermediate between zero and 
one. The critical region will thus have a finite level compressibility different from both metal 
and insulator regimes and therefore multifractal nature of the eigenvectors . 

The statistical measures for the Anderson transition in presence of a time-reversal sym- 
metry can similarly be obtained by using their equivalence to those for the Poisson — > 
GOE transition. Although much information is not available about the latter (due to the 
technical difficulties in integrating over associated orthogonal space of eigenvectors) how- 
ever some approximate results are known. For example, the R 2 for small r can be ob- 
tained by solving the eq.(6) for (3 = 1, R 2 (r, A) w (2r - l) 1/2 J 1/3 ((2r - l) 3 / 2 /3A) e r / 2A with 
J(z) as the Bessel function. For r << y/X, R 2 oc r implying level- repulsion of the same 
degree as for the delocalized regime; this is consistent with description of A as g 2 . For 
large-r, R 2 can be obtained from eq.(7) by taking R 2 (r;0) = 1 (Poisson limit), (3 = 1 and 
Mr, oo) = l-^-(r d^) (GOE limit); R 2 (r, A) » R 2 (r, oo) + r2+4 4 ^ A2 

The exact formulation for the level compressibility x f° r this case can not be obtained 
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due to lack of the knowledge of B,2(r,A) for entire energy-range. However the presence of 
a non-zero, non-unity compressibility can also be seen from the nearest-neighbour spacing 
distribution P(s) for large s values. This is because the compressibility of the spectrum is 
basically a reflection of the correlations of levels at large distances. For WD case, P(s) oc 
se -7rs 2 /4 decays as a Gaussian near the tail indicating no compressibility of levels. As for 
Poisson limit P{s) oc e~ s , the unity coefficient of s in the exponent implies a very high 
degree of compressibility. An exponential tail of P(s) with a non-unity coefficient will 
thus correspond to a a compressibility less than Poisson case; if P(s) oc e _fcs for large s, 
X ~ l/2k. The P(s) for MI transition can be given by using the one for the Poisson — > 
GOE transition @, P{s,A) = s J^ e -u 2 Ws 2 /^ 2 J Q °° dCe-« 2 - 2 « A / (sC«(A)/A) with u{A) = 



y/itU{— 1/2, 0, A 2 ), U being Kummer function and Jo Bessel function. Although this result 
is rigorous for 2x2 matrix space but is proved reliable for systems with many levels. 
(The dominance of nearest-neighbour interaction is also supported by the mapping of the 
eigenvalues dynamics to the particle dynamics of calogero Hamiltonian with inverse-square 
interaction). As can be checked, for s >> A, P(s,A) oc e~ ks with k = \/irA. Near the 
critical point, therefore, k = V7rA* > 1 and x ~ 1/2V 7rA*. For small separations, that 
is, s << A, P(s,A*) has a linear dependence on s, indicating a WD type behaviour. The 
spacing distribution for CSS is thus a hybrid of Poisson and WD type distributions; this 
result is in agreement with numerical study of 0. 

The case of the MI transition due to off-diagonal Gaussian disorder can similarly be 
studied by taking finite a^i values for all k, /; the Y for any intermediate stage of the 
transition brought about by a most general type of diagonal and off-diagonal disorder is 
given by eq.(3). As obvious from the above, the presence or absence of the off-diagonal 
disorder will only affect the form of Y, leaving, therefore, the qualitative behaviour of various 
correlators unaffected. The presence of different boundary conditions (e.g. periodic) or 
different topology for disordered lattice will also result in different F-values (as various y^s 
will be affected), leading to different level-statistics even if degree of disorder, hopping rate 
and dimensionality is same; this is in agreement with numerical observations ||13|| . 
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It should be noticed here that the expressions of x(A) as well as P(s, A) are also valid 
for Rosenzweig-Porter (RP) ensemble (< HZ >= V\ for all i, < Hf- >= v 2 for all i 7^ j and 
V\,V2 different in general), thus implying a non-zero compressibility for finite A- values in RP 
model too contrary to earlier claims |jTT|. Our result is in agreement with Dyson's conjecture 
for Brownian ensembles J7J according to which the parameter A governs a smooth transition 
in fluctuation properties; the large- s behaviour too, of P(s, A), is therefore expected to vary 
smoothly from e~ s (Poisson) to e _/3s2 (Wigner). 

To summarize, our analytical study shows that the statistical transition of levels of a 
disordered system due to change of the parameters e.g. disorder, and system size can be 
reduced to the transition with respect to a single scaling parameter. Our results also confirm 
the presence of just three type of universal statitsics for inifinite system, namely Poisson 
above critical value of the disorder, WD below critical value of disorder and a new type of 
statistics for the critical region. The critical region statistics is a hybrid of the Poisson and 
WD and shows a finite level compressibility. The universality of CSS is restricted in the sense 
that, although system-independent, it depends on the critical value A = A* too besides the 
underlying symmetry of the system. We have discussed here only the 2 nd order level-statistics 
for Anderson transition. The information about the higher order correlations can be obtained 
from the hierarchical equation (5) for the correlators or by analyzing particle correlators of 
CSM Hamiltonian. Further the mapping of the level-dynamics during Anderson transition 
to particle-dynamics of CSM Hamiltonian once more indicates the dominance of two body 
interaction in nature. 
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